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Motivated by the recent proposed models of the information engine [D. Mandal and C. Jarzynski, 

Proc. Natl. Acad. Sci. 109, 11641 (2012)] and the information refrigerator [D. Mandal, H. T. 

Quan, and C. Jarzynski, Phys. Rev. Lett. Ill, 030602 (2013)], we propose a minimal model of 
the information pump and the information eraser based on enzyme kinetics. This device can either 
pump molecules against the chemical potential gradient by consuming the information encoded in 
the bit stream or (partially) erase the information encoded in the bit stream by consuming the Gibbs 
free energy. The dynamics of this model is solved exactly, and the “phase diagram” of the operation 
regimes is determined. The efficiency and the power of the information machine is analyzed. The 
validity of the second law of thermodynamics within our model is clarihed. Our model offers a 
simple paradigm for the investigating of the thermodynamics of information processing involving 
the chemical potential in small systems. 

PACS numbers: 05.70.Ln,05.40.-a, 82.30.-b, 82.39.Fk 


I. INTRODUCTION 

In 1871 [1], in a thought experiment, James C. 
Maxwell conceived an intelligent creature, now known 
as Maxwell’s demon [2] (MD), to elaborate the statisti¬ 
cal nature of the second law of thermodynamics. In this 
thought experiment it was assumed that this demon can 
separate the fast-moving molecules from the slow ones in 
the gas contained in a box, thus apparently “violate” the 
Clausius’ formulation of the second law. Since its concep¬ 
tion, the consequence of the MD has been under debate 
by a lot of researchers, including many renowned physi¬ 
cists for many decades [2], and the MD has become one 
of the most famous thought experiments in the history 
of physics. Many interpretations on the conceptual levels 
have been put forward to guarantee that the MD will not 
lead to the violation of the second law. Nevertheless, it 
was not until 1982 when C. H. Bennett gave an expla¬ 
nation [3] based on a principle proposed by R. Landauer 
[4] (and independently by O. Penrose [5]), a consensus is 
finally reached. According to Bennett [3], information, 
or Shannon’s entropy, of the memory register of the de¬ 
mon plays a key role in preventing the second law from 
being violated in the Maxwell’s thought experiment. The 
memory register of the demon has to be erased in a cy¬ 
cle, consuming /cBTln2 energy per bit (Here is the 
Boltzmann’s constant, and T is the temperature of the 
reservoir). When the erasure of the memory register is 
included into the cycle, the total entropy production can 
never be negative. Thus the second law of thermody¬ 
namics will not be violated by the MD [3]. 

In principle, the low-entropy state of the memory regis¬ 
ter can help rectify the thermal fluctuations, and convert 
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one form of disordered energy (e.g., the heat) into an¬ 
other form of ordered energy (e.g., the work). The cost 
is that the Shannon’s entropy of the memory register will 
increase by at least the amount of the entropy decrease 
in the non-information-bearing degrees of freedom (e.g. 
the heat reservoir). In practice, however, the implemen¬ 
tation of the Maxwell’s thought experiment in a labo¬ 
ratory and the realization of the “tradeoff’ between the 
energy and the information has not been fully demon¬ 
strated in the laboratory. It has attracted increased at¬ 
tention in recent years to bring Maxwell’s thought exper¬ 
iment into practice in either artificial or natural systems. 
The candidates include systems borrowed from the quan¬ 
tum information processing [6-8] , the signal transduction 
in biological network [9, 10], the AMO systems [11-14], 
and the artificial nano-scale devices with feedback con¬ 
trol [15-30]. In almost all these models, an intelligent 
being or an external controller is involved (one exception 
is [14]). It is less interesting to have an intelligent being or 
an external controller involved when realizing Maxwell’s 
thought experiment since there is no demon in the world. 
Thus, it has become a goal for many physicists to build 
an “autonomous” information machine, which can rec¬ 
tify the thermal fluctuations and mimic what Maxwell 
has conceived but without involving any intelligent be¬ 
ing. Along this line, an autonomous information engine 
model [31, 32] and an autonomous information refrigera¬ 
tor model [33] have been proposed. Also, an optical MD, 
which can reduce the entropy of the atoms at the expense 
of increasing the entropy of scattered photons has been 
explored [11-14]. Nevertheless it is difficult to charac¬ 
terize the entropy increase in the scattered photons and 
study the second law as well as the efficiency and the 
power of the optical MD. 

In this paper we propose a minimal solvable model 
of an autonomous information pump and eraser that 
connects the information and the chemical potential: it 
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pumps molecules against the chemical potential gradi¬ 
ent by simultaneously writing information to the mem¬ 
ory register, without the expenditure of work or free en¬ 
ergy, and without the involvement of any intelligent be¬ 
ing. Namely, it is an autonomous information machine. 
We notice that there are realistic molecule pumps in bi¬ 
ological systems, such as ATPase [34-36] or ion chan¬ 
nels [37]. Usually one kind of molecules (such as Na'^) 
are pumped against the chemical potential gradient at 
the cost of dissipating the free energy of other molecules 
(such as ATP). However, in our model, instead of con¬ 
suming the chemical potential difference of another kind 
of molecules, we consume the information. The model 
is constructed based on classical enzyme reactions with 
three states [38-40], coupled to a memory register, which 
carries bits of either 1 or 0 in a stream (Fig. 1(b)). The 
enzyme undergoes stochastic evolutions, which catalyses 
the reactions and overrides the states of the bits. In cer¬ 
tain regions of the parameter space, the information ma¬ 
chine evolves to a periodic steady state and continuously 
pumps molecules, while in other regions of the parame¬ 
ter space, it (partially) erases information encoded in the 
memory register at the cost of dissipating free energy. 

This paper is organized as follows. In section II we 
introduce the three-state enzymatic system and the bit 
stream that couples to the three-state system. In sec¬ 
tion III we solve the dynamics exactly, and discuss the 
the efficiency and the power of our model. The validity 
of the second law of thermodynamics within our model 
is clarified. In section IV, we discuss and conclude our 
paper. 


II. THE MODEL 


The original enzymatic model which describes how an 
enzyme E catalyzes the reaction [38-42] is shown in Fig. 
1(a). A E molecule binds with a S molecule and forms 
a complex ES molecule. A ES molecule undergoes a cat¬ 
alytic process and becomes a new complex EP molecule. 
The EP molecule then releases the product - a. P molecule 
and the enzyme E is set free and goes into the next cy¬ 
cle. All these reactions take place under isothermal and 
isobaric conditions, so the Gibbs free energy (chemical 
potential) can be used to determine the direction of the 
reaction. The stochastic dynamics is governed by the 
master equation 

^=Koq«). (1) 

where 
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and q(t) = [qEs{t)^QE{t),QEp{t)]'^ describes the prob¬ 
ability distribution on the three states ES, E and 
EP. When the chemical potential difference between P 
molecules and S molecules 
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is fixed at a nonzero value, in the long time limit, the 
system reaches a nonequilibrium steady state. Here 
/3 = l/ksT. = {kik2kz)/{k-ik_2k-3)-, h {i = 
1, 2, 3, —1, —2, —3) describes the reaction rate constant. 
[P] and [S'] depict the concentration of P molecules and 
S molecules. In this case, the Gibbs free energy is con¬ 
tinuously dissipated into heat and released to the heat 
reservoir. The difference of the standard Gibbs free en¬ 
ergy of formation AG° of the two states connected by the 
jth reaction is determined by the dynamic parameters: 

Usually the chemical potential of S molecules should be 
higher than that of P molecules so that the substrate 
can be converted into the product. The chemical poten¬ 
tial difference is the driving force of the transformation 
from S molecules to P molecules. When the chemical 
potential difference disappears (A/r = 0), the system ap¬ 
proaches an equilibrium state. That means the forward 
and the backward reactions are balanced, and there is 
no conversion from S molecules to P molecules, or vice 
versa. 

We now combine this enzymatic model with the in¬ 
formation encoded in the incoming bit stream (Fig. 1 
(b)) and introduce the competition between the chem¬ 
ical potential difference A^ and the information (Fig. 
1(c)). Our model is very similar to the information en¬ 
gine model proposed in Ref. [31] and the information re¬ 
frigerator model proposed in [33] . The information is en¬ 
coded in the bit stream that passes through and interacts 
with the three-state enzyme system. A single bit in the 
stream, labelled with 0 or 1, can couple to the three-state 
system, and forms a composite system of six states, EO, 
ESO, EPO, El, ESI and EPl. The information machine 
is based on the six-state system and can jump within 
EPl ^ El, El ^ ESI, EPO ^ EO and EO ^ ESO 
(black solid line of Fig. 1(c)) without flipping the bit. 
This kind of evolution is simply the binding-unbinding 
process without any molecular conformational change. 
In addition, the enzyme makes a catalytic reaction from 
ES molecules to EP molecules (from EP molecules to ES 
molecules) if and only if the bit flips from 1 to 0 (from 0 to 
1), as indicated by the transformation (red dashed line in 
Fig. 1(c)) between ESI ^ EPO. Notice that there is no 
transformation between ESO and EPl, ESI and EPl, or 
between ESO and EPO. We will see that this asymmetry 
in the dynamics is essential for the pumping of molecules 
or the erasure of information. This new enzyme kinetics 
is governed by the master equation 
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denotes the probability distribution in the six states 
ESO, EO, EPO, ESI, El and EPl. All the parameters 
in Eq. (6) are the same as those in Eq. (2). The 
difference between the six-state dynamics (5) and the 
three-state dynamics (1) is that the six-state system will 
reach a periodic steady state, instead of a steady state, 
because of the interruptions of the incoming bits and the 
finiteness of the interaction time with every bit. As the 
bit stream passes by, the three-state system is assumed 
to interact with only one bit (the nearest bit) at a time. 
The duration of the interaction with each bit is set to 
be r, which is independent of other parameters and can 
be tuned. We assume that the incoming bit contains a 
mixture of I’s and O’s with probabilities pi and po- Let 

S=po-pi (7) 

denote the excess of O’s in the incoming bit stream. The 
incoming bits are statistically independent of one an¬ 
other. When a fresh bit comes in, it evolves with the 
three-state system for a time period of r, jumping among 
the six states. After each interaction period, the bit in 
the stream can be overridden due to the dynamic evo¬ 
lution (5) of the composite system. At the end of the 
period, the final state is decomposed into the bit and 
the three-state system. The final state of the bit is pre¬ 
served in the outgoing bit stream, and the final state of 
the three-state system is recycled as the initial state for 
the next period. For example, if the state at the end of 
the nth period is ESI, and the (n -|- l)th incoming bit 
is 0, then the initial state of the six-state system in the 
(n 4- l)th period is ESO. This effective transition does 
not imply an actual change in the state of a given bit, 
but simply reflects the replacement of an outgoing bit in 
state 1 with an incoming bit in state 0 [31, 33]. Given all 
the dynamic parameters k±i and S, after many periods, 
the three-state system will reach a periodic steady state. 
Finally, let 6„ and b'^ denote the incoming and the out¬ 
going state of the nth bit in the stream. The change of 
the bit during the nth interaction period is 

Axn = K-K- (8) 

Notice that Ax„ also indicates the net number of P 
molecules being transformed to S molecules during the 
nth interaction period. 

Before going to the exact solutions of the dynamics, 
we give an intuitive explanation of how the information 


tial reservoir {P constituents in the solution) to the high 
chemical potential reservoir (5 constituents in the so¬ 
lution). For this purpose, we assume all the bits in 
the incoming stream are in state 0, so that d = 1. In 
this circumstance, in every period the information ma¬ 
chine starts from one of three states EO, EPO or ESO. 
The composite system of the three-state system and the 
bit evolves for a time period of r, according to the dy¬ 
namics described by the master equation (5). The six- 
state system might jump forward and backward between 
ESI and EPO repeatedly. At the end of the interac¬ 
tion period, if the six-state system ends up in state EO, 
EPO or ESO, it must be the case that every transition 
ESI —>■ EPO is balanced by a transition EPO —>■ ESI] 
hence = 0 (8) and there is no net transformation 
between P molecules and S molecules. If the final state 
is El, EPl or ESI, we know that there is a net transition 
from EPO to ESI] then Axn = 1 (8), and the result is 
that one P molecule is transformed to one S molecule. 
This net transformation will be recorded in the outgoing 
bit stream, that a bit initially in state 0 will be finally 
rewritten to 1. In the long run, 

[t/r] 

x{t) = (9) 

n—1 

grows with time and we will get a net transformation 
from P molecules to S molecules, which is proportional 
to the I’s appearance in the outgoing bit stream. 

During this process, Gibbs free energy is generated at the 
cost of increasing the entropy of the bit stream, instead 
of the expenditure of work or free energy. Please notice 
that the mode of the energy conversion in the information 
pump is different from that of the information engine or 
the information refrigerator. In the information engine 
(refrigerator), work is generated (heat is pumped) at the 
cost of increasing the entropy of the bit stream. 

A more general case is when the incoming bit stream 
is a mixture of O’s and I’s. The excess of 0 over 1 (d > 
0), or excess of 1 over 0 (<5 < 0), resists or enhances 
the normal flux stream from S molecules to P molecules. 
The fraction of O’s in the incoming bit stream will affect 
the flow of molecules being catalysed, as well as the bit 
stream that passes by. In the following we will give a 
quantitative analysis on these two effects. 
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FIG. 1: Schematic figure of enzyme kinetics (a), the bit stream (b) and the information pump or information eraser based on 
the enzyme kinetics (c). 


III. SYSTEM AS INFORMATION PUMP OR 
INFORMATION ERASER 

A. solution to the model 

Similar to the information engine model [31] and the 
information refrigerator model [33], the evolution of the 
composite six-state system in every period can be sepa¬ 
rated into two stages: (i) the dynamic evolution under 
the governing of Eq. (5), and (ii) the projection process 
at the end of every period, which eliminates the correla¬ 
tions between the three-state system and the bit. For a 
finite period r, the three-state system will reach a peri¬ 
odic steady state after many periods (this result is guar¬ 
anteed by the Perron-Frobenius theorem [31]). In order 
to obtain the result of the periodic steady state, we fol¬ 
low the same procedure as that in Refs. [31, 33]. When 
the three-state system reaches the periodic steady state, 
we use Pq and p'l to denote the outgoing bits’ fraction 
for O’s and I’s. Let 6' = Pq — p'l be the excess of O’s in 
outgoing bit stream. The average production of I’s per 
interaction period is characterized by 

^ = (Ay) =p'i-Pi =(10) 

As described before, for each flipping of 0 —>■ 1 there 
must be an accompanied molecular conformational trans¬ 
formation of the three-state system from EP to ES, 
and accordingly a transformation from P molecules to 
S molecules. The average Gibbs free energy gained by 
the information machine per interaction period is 

AG = -$A^. (11) 

A positive AG indicates that the molecules flow against 
the chemical potential gradient, and P molecules are con¬ 
verted to S molecules continuously. 

We use the Shannon’s information entropy 

I (<5) = -pQ log 2 Po - Pi log 2 Pi 

1-d, 1-5 14-5, l-b5 (12) 


to characterize the information content of every bit in the 
incoming bit stream. Information content of every bit in 
the outgoing bit stream 1(5') could be defined similarly. 
The information entropy increase of every bit is given 
by the difference of the information entropy between the 
outgoing and the incoming bits 

AJ = 7(5') - 7(5) = 7(5 - 2$) - 7(5). (13) 

For simplicity, we consider the case where the Gibbs 
free energy difference comes only from the concentration 
difference between S molecules and P molecules. That is 
we choose k±i = 1, (AG° = 0) and let [S'] = 1 -|- e, [P] = 
1 — e, e G (—1,1). Given a set of parameters (e, 5; r), we 
obtain the exact solution of the periodic steady state (see 
Appendix A) and find the expression of the current 

4>(5,e;T) = ^[l-i77(r)], (14a) 

where 

_ 2 . (1 + 8a + 4^3^) - (2 + 7a + 4V3/3)e-2" 

® 3 - (2 + a)e-2- 

(14b) 

and a = cosh(\/3T), f3 = sinh(-\/3r). Notice that $ de¬ 
termines both A7 and AG. So we obtain the entropy 
increase in the bit stream and the free energy gaining 
in every period from Eq. (14a-14b). The term 5 — e 
in Eq. (14a) demonstrates the competition between the 
chemical potential gradient and the information. Actu¬ 
ally, we have got exact solutions for a general class of 
models, of which our current model and the model in 
Ref. [31] are two special cases (see Appendix B). 

In the parameter space of 5 and e, we plot a “phase 
diagram” (for different r the boundary of the “phase” 
are different). The information machine operates as an 
information pump when the information is the predom¬ 
inant driving force 0<e<5,<l’>0 (more S molecules 
and less P molecules) or5<e<0,$<0 (more P 
molecules and less S molecules). These conditions corre¬ 
spond to the blue regions in Fig. 2. In these regions, the 
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information machine pumps the molecules from the low 
chemical potential reservoir to the high chemical poten¬ 
tial reservoir, at the expense of increasing the entropy of 
the bit stream (writing information to the bit stream). 

The information machine can also operate as an in¬ 
formation eraser when the chemical potential difference 
is the predominant driving force, highlighted in Fig. 2 
with red. When (5 > 0, namely there are more O’s 
than I’s in the incoming bit stream, the information 
will be (partially) erased when pg > pg (more 0 is writ¬ 
ten) or pg < Pi (more 1 is written), which is equiva¬ 
lent to e > (5 or e < —[3 -I- K{t)\5/[2> — K{t)\. When 
(5 < 0, the information erasing conditions will be e < <5 or 
e > -[3-1- K{t)]5/^ — K{t)]. Working in these regions 
of parameter space, the information machine consumes 
the Gibbs free energy to reduce the information entropy 
of the bit stream. Different interaction time periods r 
correspond to different periodic steady state solutions 
(14a-14b), and hence different boundaries in the phase 
diagram (see Fig. 2). When r —>• oo, the boundary of the 
“eraser phase” will eventually approach the line e = —S. 

Our model thus reflects the competition between two 
effective driving forces. The operation regime is deter¬ 
mined by which driving force is predominant. In Fig. 2, 
one axis of the phase diagram is associated with the ran¬ 
domness of the incoming bits (i5) and the other with the 
chemical potential difference (e). When the information 
machine is operated in the shaded area, either molecules 
are pumped against the chemical potential gradient, at 
the cost of simultaneously writing information to the bit 
stream; or the information encoded in the incoming bits 
is (partially) erased, at the expense of Gibbs free energy. 
The separating line between these two regimes is given 
by i5 = e, where the two driving forces are comparable. In 
the unshaded regions, the information machine consumes 
both resources and is a dud. 


B. efficiency and power as functions of the 
parameters 

Besides determining the operation regimes of the in¬ 
formation machine, we can further study the efficiency 
and the power of the information machine for an arbi¬ 
trary given set of parameters (d,e,r). AG (11) is the 
free energy gained by the information machine during 
one period. The efficiency of the information pump can 
be expressed as 


r]{S,e;T) 


AG 

\n2kBTAI" 


(15) 


which approaches unity when d —> e, but in this case the 
power is vanishingly small (see Fig. 3). This result can 
be understood intuitively. When S ^ e, the two driv¬ 
ing forces are almost equal to each other, so the process 
is very close to equilibrium. In another word, the irre¬ 
versible entropy production is vanishingly small. Hence, 
the efficiency approaches its maximum value - unity. In 



FIG. 2: (Color online) “Phase diagram” of the information 
pump and eraser in the parameter space. The numbers indi¬ 
cate different lengths of the time periods r = 0.1,2,5,oo. 


the information eraser regime, one can define the effi¬ 
ciency as 


r]'{S,e;T) 


In 2kBTAI 

AG 


(16) 


and one can easily find that yy' —>■ 1 when (5 —> e for the 
same reason (see Fig. 3). 

The power of the information pump AG /t is also il¬ 
lustrated in Fig. 3 (we do not study the power of the 
information eraser). It can be seen that the maximum 
power is reached when the driving force of the infor¬ 
mation reaches its maximum value i5 = ±1 and when 
the chemical potential difference is in the intermediate 
regime. This result can be understood intuitively: When 
the chemical potential difference A/r is very large, the av¬ 
erage number of molecules pumped in every period $ is 
small. Hence, the power is low. When Afi is very small, 
the average number of molecules pumped in every period 
is large, but the free energy gaining A/i associated with 
one particle is very small. Hence, in both cases the power 
does not reach the maximum value. 


C. the second law 

Similar to the information engine [31] and the infor¬ 
mation refrigerator [33] , one can prove that the following 
inequality is satisfied in our model 

AG < In 2A:B'rA/ (17) 

for any given set of (e,(5, r). The equality holds only 
when e = 5 (see the SI in Ref. [31]). This is also the up¬ 
per limit of the efficiency for the information machine to 
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FIG. 3: Efficiency 77 and power t^Gjr of the information 
pump and information eraser for t — 5 and 13 = 1. 


pump molecules against the chemical potential gradient. 
When the information machine operates as an eraser, 
I{S') < I{S), the inequality is lAG] > ln2fcBr|A/|. This 
is also the upper limit of the efficiency to erase informa¬ 
tion. Both inequalities imply the Landauer’s principle [4]. 
The entropy decrease of the reservoir, including S and P 
constituents as well as the heat bath, due to molecules 
transformation can be defined as AG = —AG/\n2kBT, 
and it must not exceed the entropy increase in the bit 
stream: 


Air + A/ > 0. (18) 

If we identify the information entropy of the bit steam 
with the thermodynamic entropy, this inequality (18) in¬ 
dicates that the second law of thermodynamics is con¬ 
served in our autonomous information machine [31, 33], 
and the Landauer’s principle is equivalent to other for¬ 
mulations of the second law. 


IV. DISCUSSIONS AND SUMMARY 

In this article we propose a model of information pump 
and information eraser. This model provides a simple 
paradigm for the investigation of the thermodynamics of 
information processing based on enzyme kinetics. 


A. difference between the nsnal thermal engines 
and stochastic thermal machines 

The original development of thermodynamics found its 
inspiration in the steam engine. Nowadays the field is 
driven by the tiny molecular engines within living cells. 
Though of vastly differing scales, these engines share a 
common function: they transform energy into motion. 
For instance, ATP molecules provide the fuel for myosin 
molecules in muscle tissue to move along actin filaments, 
pulling the muscle fibers to which they are attached. 
Other motors are powered by light, by differences in 
proton concentrations or by differences in temperature 
[43, 44]. Obviously the two kinds of engines also have dif¬ 
ferent features (see e.g., section 9.6 of [45]). In the usual 
thermodynamic machines, such as a Carnot engine based 
on the classical ideal gas, the working substance for the 
thermodynamic machine is macroscopic in size, and usu¬ 
ally the ensemble average of the thermodynamic variables 
is adequate to describe the system and the thermal fluc¬ 
tuations are usually ignored. In order to convert one form 
of energy into another form, the system must undergo 
certain thermodynamic cycles. However, in a stochastic 
thermal machine, for example in the Feynman’s ratchet 
and Pawl [46], Butticker-Landauer’s ratchet [47, 48], or 
in the information engine [31], there are no thermody¬ 
namic cycles. They usually operate in (periodic) steady 
states. In addition, the thermal fluctuations are essen¬ 
tial for the operation of the stochastic thermal machines, 
which puts a constraint on the upper limit of the size of 
the stochastic machines - it must be microscopic in size 
so that it can be driven by the thermal fluctuations. 


B. competing driving forces and stochastic thermal 
machines 

Generally speaking, a stochastic thermal machine is 
a special kind of system that rectifies thermal fluctua¬ 
tions and converts one form of disordered energy into 
another form of ordered energy, by consuming certain 
kind of ordered resource, such as the information or the 
free energy. Two resources, or driving forces, compete 
with each other, generating one continuously in the cost 
of consuming another one. The driving forces can be 
the information, the temperature difference, the mechan¬ 
ical work, and the chemical potential difference. Fig. 4 
summarizes several known stochastic thermal machines. 
Besides the three information machines we mentioned in 
table I, there are three kinds of other stochastic ther- 
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FIG. 4: Stochastic thermal machines and four kinds of ther¬ 
modynamic driving forces; mechanical work W, chemical po¬ 
tential difference Afi, temperature difference AT, and infor¬ 
mation AI. Different stochastic thermal machines can be 
constructed between any pair of these thermodynamic driv¬ 
ing forces. 


mal machines: (1) the molecular motor which converts 
chemical potential difference, such as Gibbs free energy 
from ATP hydrolysis [39, 49-51], into mechanical work, 
in the form of pulling loads along microtubes; (2) Feyn¬ 
man’s rachet and pawl [46] in which temperature gradient 
competes with the mechanical work; (3) thermoelectric 
devices [52] in which the chemical potential difference 
(voltage) and the temperature difference compete with 
each other ( Seeback effect and Peltier effect). 


C. information machines, feedback control and 
Maxwell’s demon 

Among the stochastic thermal machines (Fig. 4), the 
information pump (current study), the information en¬ 
gine [31] and the information refrigerator [33] form a 
category of information machines (see Table I). They 
rectify thermal fluctuations and convert the disordered 
energy into other forms of ordered energy, such as the 
mechanical work by simultaneously consuming the infor¬ 
mation. It is worth mentioning that, they are different 
from a class of models with feedback control [18] where 
the mutual information (or correlations) between the sys¬ 
tem and the measurement apparatus is identified as the 
information content. The mutual information can be har¬ 
nessed to extract the ordered energy by rectifying the 
thermal fluctuations. In the information machine mod¬ 
els, however, the correlations between the three-state sys¬ 
tem and the bits are eliminated artificially at the end of 
every period. Also, the information machines differ from 
the usual MD-assisted thermodynamic machines in the 
following aspects [8]: In the former cases it is impossi¬ 
ble to identify the measurement process and the feedback 
control process. All processes are autonomous and simul¬ 
taneous. Thus, the intelligent being is unnecessary. In 
the latter cases, however, there must be an intelligent be¬ 
ing (or an external controller), so that the actions of the 
MD, such as the measurement and the feedback control, 
can be carried out. 


D. molecular pumps in biological systems 

In the past two decades, the thermodynamic machin¬ 
ery of life [45] has attracted a lot of attention, and studies 
in this field bloomed in recent years [38, 53-62]. Actually, 
in biological systems, instead of the temperature gradi¬ 
ent, the chemical potential is a prevailing driving force 
[63]. Living systems consume the Gibbs free energy to 
perform physical and chemical functions [45, 57, 64]. The 
chemical potential difference comes mainly from ATP hy¬ 
drolysis, or ion transport. 

In our model, the enzyme is a stochastic machine and 
the enzyme reactions are naturally built on the molecu¬ 
lar level. The enzyme system [38] is less artificial than 
the mechanical models [31-33] and has been studied ex¬ 
tensively. The experimental techniques accumulated in 
the field of enzyme kinetics will facilitate the realization 
of the Maxwell’s thought experiment. In fact, there are 
many cycle reactions that have net fluxes in cells, such as 
the phosphorylation-dephosphorylation cycle [40] driven 
by ATP hydrolysis. As long as we can design reactions 
with asymmetric kinetics like the one designed in our 
model, it is possible to build an autonomous informa¬ 
tion pump that produces molecules with higher chemical 
potential, without consuming other free energies. 


E. concluding remarks 

In summary, in this article, we propose an autonomous 
information machine that can pump molecules against 
the chemical potential gradient without the expenditure 
of work or free energy, or (partially) erase information 
encoded in the bit stream. Information is a new driv¬ 
ing force that can be used to pump molecules against 
the chemical potential gradient. We solve the dynamics 
exactly and obtain the “phase diagram” in the parame¬ 
ter space. In certain regions of the parameter space the 
information machine operates as an information pump, 
while in some other regions of the parameter space it 
operates as an information eraser. In both regions, the 
Landauer’s principle (one formulation of the second law 
of thermodynamics) is conserved. Our model offers a 
simple paradigm for the studies of thermodynamics of in¬ 
formation processing by small systems. Our investigation 
not only brings important insights to the understanding 
about the intrinsic relations between information, free en¬ 
ergy, entropy, and work, but also sheds new light on the 
ultimate realization of the Maxwell’s thought experiment 
in the laboratory. 
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Information Machines 

Input 

Output 

Information Engine [31] 
Information Refrigerator [33] 
Information Pump 

Information 

Information 

Information 

Mechanical work 

Heat flow against temperature gradient 
Molecule flux against chemical potential gradient 


Foundation of China under grant 11375012, and The Re¬ 
cruitment Program of Global Youth Experts of China. 


Appendix A: solution to the information pump 

From Eq. (6) and the assumptions in section 111. A, 
the transition matrix can be written as: 


To + Pi ■ After a straightforward calculation we obtain: 


, / F + G + SH M-2SL F - G + SH 
r=— M M + M 

\ F-G-5H M + 2SL F + G - SH 


-f 


e 

12 


F-G- H M + 2L F-G- H\ 

0 0 0 , 

G + H -F -M -2L G + H -F j 

(A4) 


/ -1 1 -he 0 0 0 

1-210 0 
0 1-e -2 1 0 

0 0 1 -2 1-he 

0 0 0 1 -2 

Vo 0 00 1-e 


0 
0 
0 
0 
1 

- 1 / 


(Al) 


with cr = e"^, F = 4 + 2a^,G = 4a^ + a‘^ + a'^,H = 
V3(cr°-cr‘^),L = 2cr2_cr‘=-cr'', and M = 4-4cr^. Based 
on Eq. (A4) we obtain: 


qP"" = 


1 

3 


1 + F 
1 

1-E 


(H-L)(S-e) 

6-G 


-he. 


(A5) 


We find that the matrix 7Z can be decomposed into 
TZ = UN~^AV. Here A = diag(0, —c, —1, —2, —3, —d) 
and N = diag(6,12c, 4, 6,12,12d), which are exactly the 
same as those in Ref. [31] while U and V are different: 


1 

X 

1 

X 

xy ^ 

X 

\ 

y 

—a 

0 

-1 

-22/-1 

-h 

1 

a-\-y 

-1 

— X 

1 

b + y 


1 

—a — X 

-1 

y 

1 

—b — X 


x~^ 

a 

0 

1 

1 

to 

1 

b 


V x-^y 

-y 

1 

-y 

X 

-y 

) 


(A2a) 


/ xy 

xy 

xy 

xy 

xy 

xy \ 

y 

-ay 

a + y 

— {a + x) 

ax 

— X 

y^ 

0 

-xy 

-xy 

0 

X^ 

y 

-y 

— X 

y 

X 

— X 

xy 

-2xy 

xy 

xy 

-2xy 

xy 

V y 

-by 

h + y 

-{b + x) 

bx 

-X / 


(A2b) 

Here a = 1 — -x/S, 6 = 1 -h v^, x = l-he, y = l — e, c = 
2 - 73, d = 2 -h v^. 

Eollowing the same procedure as that in Ref. [31] we 
find the periodic steady state by solving 

rqP®®=qP®®, (A3) 

where T = Voe^^M = VoUN-^e^^VM, and Vd = 
{1,1), M = {poI,piX)'^, with I to be 3 X 3 identical 
matrix, pq and pi are determined hy S = po — pi and 1 = 


Now we can calculate $ defined by Eq. (10). Since 
^ = p'l — Pi, the only remaining problem is to solve Pi . 
Actually we have [31] 


5 ) =VBe^^Mcf^%VB = 
Pi 


A direct calculation gives: 


1110 0 0 
0 0 0 1 1 1 




(A6) 

(A7a) 


where 


K{t) = e 


_ 2 ^ (1 -h 8a + 473/3) - (2 -h 7a -h 473/3)6-^^ 


3 — (2 F a)e 


(A7b) 


Appendix B: The solution to a general transition 
matrix 


Consider a general class of models with the transition 
matrix (6): 


= 


(-bi 
hi 
0 
0 
0 
0 


ai 0 0 

-oi -bs as 0 

bs —as — b2 02 

0 b2 —02 — bi 

0 0 6i 

0 0 0 


0 0 \ 

0 0 

0 0 

oi 0 

-Oi - bs as 
bs -03 / 

(BI) 
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We will give a general spectral decomposition of the ma¬ 
trix in form of 


TZ = UN-^AV = 


t 

t \ 

/ 



Ui •• 

U6 




i 

i ) 

V 


rig 

/ Ai 



/ <- Vl 


V 

Aey 


\ ^ V6 



where and bi {i = 1,2,3) are all real positive num¬ 
bers. The rows of V are the left eigenvectors, and the 
columns of U are the right eigenvectors, and VU = N. 
The eigenvalues are 


{0,UJi,Vi,U}2,V2,UJ3}, (B2) 


where vi ^2 are the roots of equation 

(A -|- ci)(A -|- C 3 ) — 0163 = 0, (B3) 

with Ci = ai + bi{i = 1,2,3) ; a;i_ 2,3 are the roots of 
equation 

(A -|- ci)(A -|- C 2 )(A -|- C3) — aib 3 {X + C 2 ) — a 2 bi(X + C3) 

- a3b2{X + Cl) = 0 . 

(B4) 

It can be proved that vi,V 2 ,u!i,uj 2,(^3 are all nonzero 
(0 is obviously not a root of Eq. (B3) or Eq. (B4)) real 
numbers. They differ from each other, so we don’t have to 
worry about the cases with degenerate or complex roots 
[65] . The corresponding right eigenvectors are 



= U = 


„ 

6163 “2 

0103 

"^«2 

O 1 O 3 

-"^«2 

0103 
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67^2 
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-^«2 
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0.2 

E(coi) 

-O 2 
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02 
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^2 
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— 62 

-G(c 02 ) 

^2 
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- 61(03 - 1 - UJ 2 ) 

02+61 

- 61(03 - 1 - 1 + 3 ) 

■^62 

CO 

1 

01 +ci 0 

CO 

1 

02 +C 16 

—6163 


where E(A) — (A-l-ci)(A-l- 63 ) — ai 63 ,G(A) — (A-l-ai)(A-|- 
C3) — 0163. The left eigenvectors of TZ can be obtained by 
solving TZ^'v'l = AfcV^, (fc = 1 , 2 , • • • , 6 ) in the manner of 


the conventional right eigenvector problem. The results 
are as follows: 



= V = 


+ C 3 ) - 


1 

b2 
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aias ^ 
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+ C 3 ) 
b2 
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1 

1 
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010361 
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-1 

as 
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1 
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02 
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(B5) 


r 


Then we calculate nk = v^u^: 

ni =(0103 -I- 0361 -I- bib 3 ){-^ + -^), 

O 163 0103 


n 2 i =62 [1 + 
+ 02)1 

n2j+i = 


{bi+ui)'^ F‘^{uJi) 
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010363 
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016163 


{{ci+Vj) -1-0163], 


(B 6 ) 

(B7) 

(B 8 ) 


where i = 1,2,3 and j = 1,2. The transition ma¬ 
trix in Ref. [31] is a special case of Eq. (Bl) when 
oi = 03 = 61 = 63 = 1, 02 = 1 -I- e, 62 = 1 — e. Cl = 2 
[i = 1,2,3). While the transition matrix ( 6 ) in our 
current investigation corresponds to another special case 
with with oi = 1 -I- e, 63 = 1 — e, 61 = 02 = 62 = 03 , i.e.. 
Cl = 2 -I- e. Cl = 2, Cl = 2 — e. 
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